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Fusion of two deuterons of room temperature energy is studied. The nuclei are in vacuum with 
no connection to any external source (electric or magnetic field, illumination, surrounding matter, 
traps, etc.) which may accelerate them. The energy of the two nuclei, in the system of center 
of mass, is conserved and remains small during the motion through the Coulomb barrier. The 
penetration through this barrier, which is the main obstacle for low-energy fusion, strongly depends 
on a form of the incident flux on the Coulomb center at large distances from it. In contrast to 
the usual scattering, the incident wave is not a single plane wave but the certain superposition of 
plane waves of the same energy and various directions. As a result of interference of them, the wave 
function close to the Coulomb center is determined by a cusp caustic which is probed by de Broglie 
waves. The particle flux gets away from the cusp and proceeds to the Coulomb center providing a 
not small probability of fusion (cusp driven tunneling) . 

PACS numbers: 25.45.-z, 03.65.Xp, 03.65.Sq 



I. INTRODUCTION 

The aspects of nuclear fusion can be found, for in- 
stance, in Refs. fJ-TTol] and references therein. Here 
we outline principal phenomena associated with nu- 
clear fusion. The main difficulty is getting the nuclei 
close enough to fuse since they should overcome a high 
Coulomb barrier. 

There are two ways to pass the Coulomb barrier, to 
accelerate the nuclei up to a high energy comparable with 
the barrier height (of the order of 1 MeV) or to pass the 
barrier via quantum tunneling. When the energy is not 
high, the probability of tunneling of the nuclei through 
the Coulomb barrier is extremely small according to the 
theory of Wentzel, Kramers, and Brillouin (WKB) [Tl| . 
So only high energy nuclei can fuse. 

Therefore, high energy nuclei is the leading idea of fu- 
sion technique. We mention an acceleration of initially 
cold deuterons by a strong electric field using a pyroelec- 
tric crystal 

It is supprizing to claim that two bare deuterons, that 
is isolated from everything, of room temperature energy 
are able, in principle, to penetrate the Coulomb barrier 
with a not small probability and to subsequently fuse. 
This statement is counterintuitive and against common 
sense which is based on the usual underbarrier physics. 

We really deal with two low-energy nuclei in vacuum. 
There are no external sources (electric or magnetic field, 
illumination, surrounding matter, traps, etc.) which may 
accelerate them. The energy of the two nuclei, in the 
system of center of mass, is conserved and remains small 
during the motion through the Coulomb barrier. This 
is a substantial difference from usual schems to push in 
action a mechanism of low-energy fusion by some local 
heating or acceleration of nuclei. 

The conventional scattering problem is a study of re- 
flection of the incident flux which is a plane wave in 
the z direction coming from large distances [111 ]. In the 
Coulomb field scattering relates to the Rutherford for- 
mula and the wave function is exponentially small at the 



center. 

What happens when the incident flux at large distances 
is not just a single plane wave? 

If the incident flux is a superposition of plane waves 
with the same energy and of various directions around 
the z axis, this can substantially modify the wave func- 
tion close to the Coulomb center because of an interfer- 
ence of the partial waves. Usually in that region the wave 
function exponentially decays along the z axis which is 
the tunneling direction. Due to the interference the wave 
function can be localized in the vicinity of the z axis 
where it weakly decays toward the center. According to 
general arguments, the kinetic energy k 2 + k 2 (for sim- 
plicity we consider the axially symmetric case), where 
r = \J x 1 + y 2 , is negative under the potential barrier. 
The decay in the tunneling direction is determined by 
the imaginary part k z = i\k z \ (the less \k z \ the more tun- 
neling probability). This means that the enhancement of 
the imaginary part | k r | , produced by the localization near 
the z axis, may reduces \k z \ down to zero and |i/>(0,z)| 
may not exponentially decay toward the Coulomb center. 
This is a redistribution of the parts in the total kinetc en- 
ergy. 

Such a state, when the wave function decays away from 
the z axis, is not physical due to the logarithmic singular- 
ity of this solution at small r [ll[ . This solution should 
be supplemented by the second branch to cancell the sin- 
gularity. But the second branch, with the opposite sign 
of k r , exponentially increases away from the z axis and 
this destroys the scheme. 

However when there is a continuous set of wave func- 
tions with various k r it may happen that branches of the 
second type cancel each other due to interference outside 
a close (non-semiclassical) vicinity of the line r = 0. In 
the semiclassical language, that set of functions does not 
correspond to a saddle point in the integration on a pa- 
rameter which marks different functions. So the state, 
localized near the z axis and weakly decayed along it (a 
channel), can be realized when the incident flux is not a 
single plane wave but is a continuous set of waves of the 
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same energy. A motion along the channel provides a not 
small fusion probability. 

One can also propose, in a less "scientific" manner, a 
general qualitative explanation of the phenomenon based 
on global properties of particle flux. It is given in Sec. Ixl 

In this paper the semiclassical method, based on 
Hamilton- Jacobi equation [TJ [T^], was used. The cru- 
cial role is played by caustics (surfaces of reflection of 
de Broglie waves of a singular shape in the form of 
a cusp which gives rise to the above channel. The flux 
gets away from the cusp to the Coulomb center along the 
channel. 

According to common sense, a tennis ball cannot pass 
through a brick wall. This correlates with the usual un- 
derbarrier physics led by the philosophy of addition of 
probabilities but not amplitudes. The latter may result 
in an unexpected scenario of barrier penetration when 
a particle passes across a classical barrier. This phe- 
nomenon was already proposed in Refs. [3 [HI]- The 
present paper develops the ideas of interference in tun- 
neling. 

Multi-dimensional tunneling was studied, in particu- 
lar, in Refs. [l6l - l24j |. Influence of interference on tun- 
neling th roug h nonstationary barriers was studied in 
Refs. [25j, |26[. The role of interference in a decay was 
studied in Ref. [I?}. Enhancement of tunneling across 
static not one-dimensional barriers due to interference 
was pointed out in Refs. [bH. liH [28jj . 

Experimental conditions to form the required incident 
flux are briefly discussed in Sec. [Xj 
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is supposed to be large which corresponds to semiclassical 
approximation. In this case one can ignore the last term 
in Eq. ([5]) and it goes over into the equation of Hamilton- 
Jacobi when a is the classical action divided by HB [ill; 

m. 

In a spherically symmetric case the first term in Eq. ([5]) 
is just (da/dR) 2 /4. The wave function exponentially de- 
cays inside the Coulomb barrier. Accordingly, the imag- 
inary part of a, related to a motion under the barrier, 
is 
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where we ignore the distance i? which is short compared 
to scales of the Coulomb motion. Eq. ([7]) follows from the 
WKB approximation The probability of fusion is one to 
overcome the Coulomb barrier 
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II. GENERAL VIEW 

Two examples of nuclear fusion are below, 
deuterium-deuterium fusion occurs by the scheme 



lu- 



lu 



3.27 Mev. 



The 



(1) 



The deuterium-tritium fusion releases more energy 

\H + \H -> \He + n + 17.59 Mev. (2) 

The problem can be separated by two steps (i) overcom- 
ing the Coulomb barrier by two nuclei with masses Mi 
and M2 and (ii) a subsequent nuclear fusion at the short 
distance Rq ~ 10~ 13 cm. The former is a motion of a 
particle with the reduced mass M — M\Mij(Mx + M 2 ) 
described by the Schrodinger equation 
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where R is the inter- nuclear distance [U| . Below we mea- 
sure length in the units of 2e 2 /E. The wave function can 
be written in the form 



tp = cxp(iBa), 



(4) 



One can estimate the fusion probability ([8]). for exam- 
ple, for the reaction (JXJ) - In this case the reduced mass 
is M — Mo/2 where the deuteron mass is defined as 
M D c 2 ~ 1.87 x 10 9 eV. Bohr radius is h 2 /Me 2 ~ 2.88 x 
10~ 12 cm. At the energy of two deuterons E — T, where 
T = 300 K ~ 2.58 x 10~ 2 eV relates to the roo m tempe r- 
ature, the de Broglie wave length is A = 2irh/ \J1MdE ~ 
1.26 A. Two deuterons enter the Coulomb barrier at the 
distance e 2 /E ~ 557 A between them. With these pa- 
rameters one can estimate ttB ~ 6174.9. This corre- 
sponds to the probability w ~ io -2682 which is in accor- 
dance with usual estimates for low-energy fusion. 



III. SCATTERING BY THE COULOMB 
CENTER 

The above estimation of nonphysically small probabil- 
ity of low-energy fusion corresponds to the usual WKB 
approximation. For a spherically symmetric incident flux 
the problem is generic with one-dimensional one since 
radial and spherical parts of the wave function are sepa- 
rated. 

There is another remarkable case of variables separa- 
tion in the Coulomb field. It occurs in parabolic coordi- 
nates 
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FIG. 1: The flux on the Coulomb center (zero coordinates) is 
reflected from the caustic surface defined by Eq. (|16l) . Four 
trajectories from the infinite set are shown. 



£ = \J r 2 + z 2 + z, ?y = \J r 2 + z 2 — z. (9) 
In this case the total wave function has the form 



(10) 



where <P is the azimuthal angle and m is the magnetic 
quantum number. The use of spherical coordinates is less 
convenient for our purposes due to the discrete summa- 
tion on azimuthal quantum number I as in the problem 
of Regge poles _29| . Below we consider an axially sym- 
metric case when m — and the Schrodinger equation 
has the form 
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where (3 is a constant connected with the variable sepa- 
ration. 

Each of Eqs. (fT2j) and (IT3l) have two independent so- 
lutions. When £ is not small (r is not small at z < 0) 
and r\ is not small (r is not small at < z) one can use 
the semiclassical approach for the wave functions which 
is not valid at very small r only [ll|. In this limit two 
branches of / have the asymptotic form 



/ T (£,/3)~exp TiB [ d^Jl-^- 
and analogously two branches of tp are 



Vt^iP) ~ ex P 
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Withing the semiclassical accuracy we do not specify pre- 
exponential factors. 



The velocity field can be studied by Newtonian trajec- 
tories. According to classical mechanics [12J, velocities in 
the £ and 77 direction are proportional to y/l — (1 + /?)/£ 
and \/l + B/rj respectively. This sets a velocity field in 
the plane {£,77}. Each point in the plane {£,77} belongs 
to one or a few classical trajectories rj(£). 

In the usual scattering problem the incident flux is 
solely a plane wave from large positive z. This situation 
corresponds to B — —i/B (f3 ~ in the semiclassical 
approximation) when the exact solution of Eq. (|13p is 
ip + = cxp(iBri) [Tlj . The combination f_(£)(p + (rj) re- 
sults in the incident plane wave ip ~ exp(— 2iBz) at large 
positive z. Analogously, the combination /+(£) <p+{i]) 
leads to the scattered wave ip ~ exp(2iBv / ?' 2 + z 2 ) far 
from the center. 

Classical trajectories are shown in Fig. [TJ They are 
reflected from the caustic surface where the velocity 
d£/dt = 0. See also [3C}. This happens at £ = 
1 + /3 ~ 1. The caustic surface in Fig. [TJis axially sym- 
metric, since m = 0, and, as follows from Eq. ([3]), is given 
by 



2z = 1 



(16) 



Along the caustic the normal momentum is zero. 
The tangent momentum is real and determined by 
sjl + B/r/ ~ 1. The wave function exponentially decays 
inside shadow region in Fig. [TJ 



IV. FORMULATION OF THE PROBLEM 

Below we explore the superposition of functions 
and (1151) with different B 



MZ,l) = J dPc I (fi)f^,p)[<p-(r,,P) + ( p+(r,,p)] 

(17) 

Mt,rj) = J d8c n {/3)f T (r},-l-/3)[<p + (Z,-l-8) 
+<p_(Z,-l-P)} (18) 

with certain weight functions cui(B) specified further. 
The total wave function is 



ip(£, 77) = ipi((„ 77) + ipn($, »?)■ 



(19) 



As shown below, the function tpj is exponentially small 
at z < and the function ipn is exponentially small at 
< z. The functions (p^irj) contain the singular part 
in 77 at small argument. The same is valid for functions 
/=f(0 PJJ- The functions </ 3 t( 7 ?) are chosen in the way 
that the non-physical singularity is cancelled in the com- 
bination <fi-(rj) + <p+(r]) which is finite at small 77. It is 
not difficult to show that at large 77 the preexponential 
factor of <p- is C /^/fj and the preexponential factor of (p+ 
is iC/y/rj where C is a constant. According to Eqs. p7|) 
and (|T5| , the main parts of the total wave function pTJ)) 
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are not singular close to r = Q at all z. Singularities of 
the exponentially small functions (tpj at z < and rpu 
at < z) are not important since they are compensated 
by exponentially small corrections of the main parts. 

The interference, described by Eqs. (IT71) and (TT%)) . of 
partial waves with different /3 strongly depends on forms 
of weight functions ciji((3). In particular, these func- 
tions define an incident flux far from the Coulomb cen- 
ter. Details of this flux, and therefore ci.n(f3), are set by 
experimental conditions. Below we use the forms 



c i,ii{0) ~ ex P (TiBa/3) . 



(20) 



where a is a fixed positive parameter. The related inci- 
dent flux is discussed in Sec. I PCI 

According to Eqs. ([14")) and flTSj) . there are four (t — 
1, 2, 3, 4) semiclassical branches of the wave function 



di exp(iBai)d/3, 



(21) 



where the preexponential factors <Zj are not important 
for the semiclassical approach used, are defined by 
the following 



gration limits only. For <ji, for example, we have 




The Cartesian derivatives 
2 



and 
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are proportional to classical velocities. 

Now one can outline the problem. From Eq. (I23|) and 
its analogs one has to determine /3(£, rj) and to insert it 
into Eq. (|22[) and its analogs. It will be a semiclassical 
solution to be found. 



CAUSTICS 




02,3 = T / \ 1 1 t~ d4i + / 4 / 1 + — e% - ia/3 

Ji+p V si io V m 

(22) 

and CTi t 4 + ia/3 = —((73,2 + ia/3). 

The large parameter B provides validity of semiclassi- 
cal approximation and the /3-integration can be done by 
the saddle method. This means that for each £ and rj 
one can determine the certain /3(^, 77) from the condition 
d<7i/d[3 = 0. For a\ this condition reads 




1+0 \i+p ' \\p ' - 1 v^J exp(a) ' 

(23) 

The analogous condition for 02 differs from (|2"3"|) by the 
sign at y/rj/p. The condition for 03 differs from (|23|) by 
the sign of a. The condition for a 4 is obtained if to make 
the both changes in Eq. (l23l) . 

The above semiclassical solutions also follow from the 
equation of Hamilton- Jacobi. Eq. ([5]) goes over into this 
equation if to drop the second derivative. The Hamilton- 
Jacobi equation in parabolic coordinates 



V 



do 
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1 



Z + V 



(24) 



allows separation of variables used above. Equations 1 
for ai with the conditions (1231) are equivalent t o a g eneral 
integral of the Hamilton- Jacobi equation (|24|) |12j |. 

To obtain derivatives, dai/d^ and dai/dr] one should 
differentiate in Eq. (f2"2"]l with respect to the upper inte- 



A condition of applicability of semiclassical approx- 
imation is well known. A wave length should be a 
smooth function of coordinates [Ijj . I n other words, 
quantum corrections to the Hamilton- Jacobi equation 
flM}, which contain (l/B)d 2 a/d£ 2 and (1/ B)d 2 a / 'd?f , 
should be small. There are various situations when those 
derivative are not small. First, it happens near a classi- 
cal turning point where classical momenta (~ da /d£) are 
proportional to square root of the distance to this point. 

In addition to that, the classical momentum can be 
not zero but nevertheless its spatial derivative turns to 
infinity. As follows from Eqs. (|2"5"|) and (|21)|) . this may 
happen when d/3/d£ or 80 /dr) become large. This con- 
dition specifies a certain surface, called caustic, in the 
three-dimensional space where classical trajectories are 
tangent to it [l3j]- In our case caustic surfaces have axial 
symmetry with respect to the z axis and one can study 
just caustic curves in the {r, z} plane. The caustic con- 
dition can be obtained from Eq. (|2"3")) by the formal con- 
dition d£/d(3 = at a fixed rj. After a little algebra the 
caustic condition reads 



2(1 + 0) cosh 2 a - sinh 2 a 



2C 

1 + /3 



± sinh ay sinh a + 4/3(1 + 0) cosh a, 

2rj 22 
— = 2/3 cosh a + sinh a 



(29) 



=F sinh ay sinh 2 a + 4/3(1 + 0) cosh 2 a , 

which determines the caustic form t/(£) if to exclude /3 in 
these equations. 



5 



The upper sign in Eq. (|29|) relates to the branch o\ and 
produces the caustics A^ and Bi in Fig. O The lower 
sign in Eq. (|29|) relates to the branch 03 and results in 
the caustics A3 and B3 in Fig. [5] 

The caustic At at small r corresponds to small j3 and 
to large j3 at large r. One can easily obtain the shape of 
the caustic At in limiting cases 



4^2 
sinh 2a 



sinh a 



3/2 



(z - z s ) < z s , (30) 



1 - 



2z s 



z tanh a 



z s < z, 



(31) 



where 2z s — cosh a. The caustic i?i corresponds to 
/3 < -1. At small r it should be (-1 - ft) < 1 and at 
large r the parameter (— /3) is large. The form of the 
caustic B\ in limiting cases is 



sinhoVl + 2z, (2z + 1) < 1, 



sinh ( 



2z s 



: tanh < 



1< z. 



(32) 



(33) 



The pair of caustics {A^^B^} is a mirror reflection of 
the pair {At, B\} with respect to the x axis as shown in 
Fig. [2] The caustics A\ and ^3 have a cusp shape close 
to the points z = ±z s . 

Along caustics Im a = const and the momentum is 
tangent to the caustic, that is, along each caustic, 



dRea\ _ / dRea\ dr 
\ dz J dz 



VI. DISTRIBUTION OF THE FLUX fa 
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FIG. 2: Distribution of the particle flux associated with the 
wave function 4>i- Caustic surfaces are axially symmetric with 
respect to the z axis. The curves represent intersection of 
caustics and the {x, z} plane. The caustic A\ touches the z 
axis at the point z s . The caustic B3 intersects the z axis at 
the point z — 1/2. On the dashed caustics fa is exponentially 
small. On the curve R the branch fa is converted into ^3. 



which are valid close to the caustic At, that is, when the 
square root in Eq. ([Bf)]) is small. Two signs of the square 
root provide two opposite velocities normal to the caustic 
At from the classical side. Close to the z axis 



r sinh a 



Ptb = ^—^, r 2 < (z - z s ) 3 < 1 



Now it easily follows from Eq. (|23|) that 



(3ta = 



(37) 



!2z sinh a + y/2z — 1 cosh ( 



(38) 



z — z s rzv2tanha 9 N o 
r — , r 2 <(z-z s ) 3 . 

Z s Z s yj Z Z s 



Pib = 



As one can see, Eqs. (|38|) go over into Eqs. (|37|) when z 
is close to z s - 



In this section we specify the velocity field associated 
with the wave function ipj. First, we focus on the region 
(z — z s ),r -C 1 where the caustic At in Fig. [2]is about to 
touch the z axis. In this case the parameter j3 is small 
and Eq. (|23p can be written in the form 



01 



3/2 



Tt +r 



— tanh a = 0, 



(35) 



which is cubic with respect to y/fit. We ascribe the index 
"1" to f3 to emphasize its connection to the branch at- 
Two physical solutions of Eq. (|3"5j) at (z — z s ) <C 1 are 



'la.b 



3z. s 



it 



3r sinh 2a 
4^ 



(36) 



A. Incident flux 

Now one can analyze what happens to the incident flux 
(from the right in Fig. [5]). It is described by the action 
o\ a where the index a means that one has to substitute 
Pta in Eq. ^ for at- As follows from Eqs. ((23]), ([271) . 
and (1281). at z not too close to z s 



data 
dz 

data 
dr 



-2\/l-_L, r< l, Zs<z 



(39) 



1 

2z 



A 



^2z sinh a + VI — 2z cosh a 
2(z - z s ) 



(40) 
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It is instructive to consider a\ on a classical trajectory 
which is just a curve in the entire space. The trajectory 
la, associated with the branch a\ a , is shown in Fig. [2] 
Trajectories of the type la do not intersect the z axis 
since, according to Eq. (|40[) , the normal to this axis mo- 
mentum becomes zero on it. 



for ip2 does not have a solution at that region. This holds 
in the semiclassical limit only, that is not too close to the 
z axis. To the left from the curve R the branch ip2 exists 
close to the z axis only where it serves to compensate the 
logarithmic singularity of tpi. So this is a realization of 
the state localized in the vicinity of the z axis mentioned 
in Sec. H 



B. Flux reflected from the caustic A\ 

The trajectory la is reflected from the caustic A\ and 
goes over after the reflection into the trajectory 16 as 
shown in Fig. [2] The related branch is crib, that is ci 
with the above ftu. As follows from Eqs. (|2li|) . (|T7|> . and 
(EH), 



da 



16 



dz 

do-ib 
dr 



-2 tanha 



I2(z - z, 



: V 2z{z — z s ) 
rZe tanha 



r < 1, 



z(z ~ z s ) 



(41) 



According to Eq. (|41l) . the trajectory 16 intersects the z 
axis since the momentum, normal to this axis, is finite. 

After intersection of the z axis, the trajectory 16 turns 
to the trajectory 2 marked in Fig. [51 The related da-ijdz 
and da2/dr differ from Eq. (j4TJ) by the sign of square 
roots. This description of the branch 2 is valid for r<l 
and z s < z. At 1/2 < z < z s and r< 1, da^jdz has the 
same form ([3"§f but dtj2 / dr differs from (|4*IH) by the sign 
of the second square root. 



C. Branch conversion 

Let us consider one of the functions ipj (TTT]) which con- 
tains f—(£)[ip—(rj) + (p+(r])]. The sum of two functions 
has no logarithmic singularity at small 77. In terms of the 
notation ([2"Tj) . this function is if>i + ip2- As follows from 
Sec. IVI1 the two branches, ipi and ip2, exist at small r 
and 1/2 < z. More precisely, they exist to the right of 
the certain curve R in Fig. [5] where two branches, ip2 and 
-03, merge due to the condition 8(12,3/ dt; = 0. 

The condition of the curve R, £ = 1 + /3, and equations 
of the type (|23l) but for i — 2, 3 yield ft — 77/ sinh 2 a. As 
a result, the curve R in Fig. [5] is determined by the form 
^r(v) = 1 + 7 l/ sinh 2 a or 



sinh a 



sinh 



-\]{1 - 2z)(sinh 2 a 



2z) 



(42) 



We consider 1 < sinh a. Close to the curve R, tangent 
derivatives are 8(72,3/ drj = 1/tanha. The derivative 
da/d^ ~ (£ — £r) corresponds to -02 to the right of the 
curve R (£r < £) and to "03 to the left of the curve R 

In contrast to "0i , the branch ip2 disappears to the left 
from R since at this region there is no saddle point for it 
in the /3-integration. In other words, the equation (|23[) 



D. Crossing of the plane {z — 0} by the flux 0/ 

As follows from Eq. (|2"3")l , on the plane {z — 0}, that is 
at £ = 77 = r, the parameter /? satisfies the equation 



1 + 2a//3(1+/3) sinh a = -/3(1 + ft) cosh 2 a (43) 



and, according to Eqs. ([251) 

<9o-i i 
— — = — tann a -= , 

dz Jr 



0, r < 1 



9(7 1 

9r 



— tanh a H -= . 

\ r 



(44) 



We see that the wave function ipi increases with z near 
the line z = 0. Close to the left cusp, {z — —z s ,r = 0}, 
the parameter ft = — 1 and, as follows from the relation 
([221 for cti, 



7T.B 



|V>j(z = -2 s ,r = 0)| ~ \ip!(z = z s ,r = 0)|exp 



The exponential comes from the 77-integration from zero 
to one. 

The caustics A\ and B3 are connected by classical 
paths, as shown in Fig. [2] For this reason is the 
same (with the exponential accuracy) on these caustics. 
The analogous statement is valid for caustics ^3 and Bi , 
indicated in Fig. [5] by the dashed curves, where \ipi\ is 
exponentially small according to Eq. ([45)) . So expo- 
nentially decays at negative z. 



VII. DISTRIBUTION OF THE FLUX ip n 

As follows from Eqs. (|4"Tj) and (|44p . the classical action 
becomes singular close to the center of the Coulomb field. 
At this region the semiclassical approximation is not true 
and the branch ipi is hybridized with 4>n- The branch 
4>n is determined by Eq. I|18p. 

We do not analyze here all details of i/jii which is sim- 
ilar to ipj. The caustics for the function ipu are shown 
in Fig. [3] where ones with exponentially small ipn are re- 
lated to dashed curves. The caustics are of the same form 
as in Fig. O The analogue of the curve R of Fig. [5] exists 
in Fig. [3] in a symmetric way at negative z (not shown 
in Fig. [3]). For the branch 16 of ipi at z s < z currents j z 
and j r are negative. For the analogous branch of ipn at 
2 < — z s the current j z is also negative but j r is positive. 



FIG. 3: Distribution of the particle flux associated with the 
wave function %j)n. On the dashed caustics ipn is exponen- 
tially small. 



This is indicated in Figs. [2] and El For the both cases the 
function /3(£, ij) is the same. The relation, analogous to 
Eq. (|45)) , for the branch ipu contains opposite sign of the 
exponential. So the branch ipu exponentially decreases 
at positive z. 



VIII. HOW THE PARTICLE APPROACHES 
THE COULOMB CENTER 

The flux on the Coulomb center comes from the right 
(-0/) and goes to the left (ipn). The functions tpi and tpu 
are hybridized in a vicinity of the Coulomb center. The 
total wave function (| 19(1 is a constant (with the exponen- 
tial accuracy) on all caustics in Figs. [5] and [3J The curves 
of classical velocities in Figs. [5] and [3] are symmetric with 
respect to the x axis. The wave function exponentially 
decays inside shadow sides of caustics. The shadow side 
of the caustics A\ and B\ is between them. The same 
is valid for the caustics A 3 and B 3 . The branches 0x2,3 
describe the whole flux distribution. The branch 04 does 
not contribute to the semiclassical wave function since a 
proper saddle point is absent. In other words, all clas- 
sical paths, described by C4, interfere down to zero not 
very close to the z axis. 

Since we are interested in the region close to the 
Coulomb center details of the incident flux, reflected from 
the caustics B\ and A3, are not important for this pur- 
pose. Therefore we focus on a\^ at small r. 

The branches a\ a and 02 decay inside the shadow of the 
caustic B3 approaching the Coulomb center. According 
to Eq. Igrj), at z s > z ~ 1 



dau 
dz 

dau 
dr 



-2 tanh a - 



ir 

z 



2z(z s - z) 



r < 1 



'2(z s — z) rz s tanha 
zz s z(z s - z) 



(46) 



We see that the branch an, also decays inside the shadow 
of the caustic A\ and dlmau/dz is reduced close to the 







FIG. 4: Curves ci, C2, and others of constant l^l are also 
directions of momenta. The semiclassical approximation is 
not valid below the dashed curve described in the text. The 
left part of the interrupted plot is related to small z. Inset 
shows in the {x, z} plane how the flux gets away from the cusp 
and proceeds to the Coulomb center (cusp driven tunneling). 



z axis (proportional to r). This is a manifestation of the 
cusp formation. 

At r ~ z <C 1 the equation (f2"5|) can be easily solved 
since in this case the parameter 1 + /3 is small. The result 
is 



daib 
dz 

daib 
dr 



2i 



tanh a 



- V 



tanh a 



z < 1 



(47) 



which also corresponds to decrease of the wave function 
inside the shadow region. The continuity equation divj = 
in semiclassical approximation is equivalent to 11] 



(VImcr) (VRecr) = 



(48) 



that is the momentum VRecr is directed along a curve of 
constant as in Fig. |U One can check that expressions 
(HI and (g7]) satisfy the condition (|48|) . 

In Fig. 2] the local angle \ between the tangent di- 
rection of the curves and the z axis is given by the 
relation 



<9Imer 



tan% 



dlma 



dr dz 
As follows from Eqs. (gB]) and (g7]), at r < 1 



(49) 



tanx 



t an \ = 



2z{z s - z) 



V r 2 + z 2 — z 



r <C z < z s 
1/2 



+ z £ + z 



r ~ z < 1. (50) 



According to these equations, the curves c\, C2, and oth- 
ers in Fig. @] are directed by 45° with respect to the r 
axis at z = 0. At z ~ 1 these curves are almost parallel 
to the z axis. 
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Another important aspect is a border of applicability 
of the semiclassical approximation used. This approxi- 
mation holds when the term with quantum corrections 
(one with second derivatives in Eq. (J5J) ) to the Hamilton- 
Jacobi equation is small. Substituting expressions (|46|) 
and (|47|) into Eq. <j5j> , one can conclude that the semiclas- 
sical approximations is violated below the dashed curve 



in Fig. @] which is 



1/B 2 atz< 1/B 2 



at 1/B 2 < z < 1/2, and r ~ 1/B at larger z. In phys- 
ical units the scale 1/B corresponds to the de Broglic 
wave length defined in Sec. HU The scale 1/B 2 is Bohr 
radius (Sec. [TT]) which is the shortest spatial scale of the 
Coulomb problem. As it should be, the semiclassical ap- 
proximation is not valid too close to the z axis (Tlj . 

It follows from Fig. 0] that the curve C2 is entirely in 
the semiclassical region and approaches to the caustic 
A\ . But the curve c\ enters under the dashed curve. 
The distribution of |0| in Fig. @] relates to the branch 
0i which has the logarithmic singularity at r = 0. To 
compensate the singularity, as discussed in Sec. EH the 
branch -02 should be added which is also singular at r = 
but tends to zero away from the z axis. 

The inset in Fig. @] shows the flux flow in the shadow 
region of the caustic A\ . According to general properties 
of caustics, \ decays at the shadow side away from A\. 
Arrowed curves close to the caustic A\ in the inset in 
Fig. prelate to constant |0i|. The closer the curve is to 
A\ the larger j^il is. The same is valid if to approach the 
z axis at the region left from the cusp, where the modulus 
of the wave function reaches its maximal value on the z 
axis. That value is kept along the whole segment of the 
z axis, from the cusp to the center. This can be called 
cusp driven tunneling. 

The caustic A\ in Fig. 2] is associated with the branch 
lb of -0/ . The analogous branch (related to 0//) is con- 
nected with the caustic A 3 and supplements Fig. [4] by the 
mirror reflection. In this way the flux is continued from 
positive to negative z. So the particle flux in Fig.(2J from 
the caustic Ai to A3, is localized at the close vicinity of 
the z axis. 

The tunneling probability can be defined as the ratio 
of densities at the center and at large distances 



w 



^(ro,0) 



V>(ri,oo) 



(51) 



The parameters ro and r\ are indicated in Fig. 01 The 
point {ro,0} is not exactly at the center but when ro ~ 
1/B 2 the expression (fSTj) does not differ from the exact 
probability in the exponential approximation. The tun- 
neling probability has the form 



exp{-2BIm [cri(r o ,0) - ai(r u 00)]} 



(52) 



The infinity point can be substituted by one on the caus- 
tic Ai as shown in Fig. [4] since to the right from the 
caustic a is real. 

If in the expression ([5"Hj) 1/B 2 < r < r^ and 1 /B < r\ 
then one can use the semiclassical approximation along 





1 26 f 



(b) 



FIG. 5: (a) Conical particle flux from the right in the absence 
of the Coulomb potential, (b) Split of the conical flux by 
the Coulomb field at large z. Trajectories are reflected from 
caustics which are continuations of ones in Fig. [2] The dashed 
parts are caustic shadows. 



the path (two thick lines) between the two dots in Fig. [4] 
For convenience, we take in addition ri C 1. Then the 
first equation (|47| simply gives 

Im[cri(ri ) oo)-(T 1 (ri,0)] = -2 v /n. (53) 
The second equation (|47l) produces 

lm[ai(r 1 ,0)-a 1 (r (h 0)}=2Vn-2V^- (5 4 ) 

By means of Eqs. (I5^|) - (fM]) the ratio of densities at 
points {ro,0} and {ri,oo} (|5ip can be written in the 
form 



w 



exp (-4B^) . 



(55) 



Equation (f55j) holds when 1 /B 2 < ro and formally that 
expression is exponentially small. It is clear that one can 
continue ro down to the border of applicability of the 
result (|55|) . namely to put ro ~ 1/B 2 . One can conclude 
from here that tunneling probability is not exponentially 
small. 



IX. INCIDENT FLUX FAR FROM THE 
COULOMB CENTER 

In the absence of the Coulomb potential the flux from 
the left to the right relates to a = — 2z and consequently 
to the plane wave = exp(— liBz) in the whole space. If 
we add the Coulomb potential that plane wave is hardly 
influenced far to the right but becomes strongly violated 
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at finite distances. In particular, the caustic is formed as 
shown in Fig.[TJ This situation corresponds to j3 = —i/B 
as mentioned in Sec. Mil 

Suppose again that the Coulomb potential is absent 
but the incident flux from the right is not a plane wave 
but one related to a = — ztanha ± rj cosh a. The cor- 
responding exact wave function is expressed through the 
Bessel function 



if) = exp {—2iBz tanh a) Jo 



/ 2Br 
\ cosh a 



(56) 



In this case the velocity distribution is shown in Fig.[5la). 

If we add the Coulomb potential the velocity field be- 
comes strongly deformed. The new velocity distribution 
is illustrated in Fig. [5jb) where 



cot 9q = sinha. 



(57) 



It is easy to qualitatively understand the features of that 
distribution. The flux, locally parallel as in Fig. 03a), is 
reflected by the Coulomb force upward when the polar 
angle 9 (tan 9 = r/z) exceeds 9 . This flux is reflected 
downward when 9 is less than 9q. There is the shadow 
region between these limits as in Fig. [5fb) . 

If n is a direction normal to a caustic surface then the 
wave function decays inside the shadow region as 



ip ~ exp 



(n/l) 



3/2 



(58) 



where / ~ (A 2 ^) 1 / 3 in physical units. Here A is the wave 
length defined in Sec. [HI Eq (|58")) is a usual form for caus- 
tics [Hj]. As follows from Eqs. (pH"]) and the distance 
between caustics A\ and B\ at large z is proportional, in 
physical units, to y/ ze 2 /E which is larger than I, So 
there is a real shadow between the caustics. 

We see that properties of the incident flux at large 
distances strongly determine the wave function at the 
Coulomb center. 



X. DISCUSSIONS 

Probability of tunneling across a one-dimensional 
static potential barrier is exponentially small when the 
barrier is almost classic. The wave function decays in- 
side the barrier since in the classically forbidden region 
the wave vector is imaginary. 

When dimensionality of a problem is higher than one 
the situation can be more complicated since the wave 
vector has more than one components in space. The sum 
of squared components, which is a kinetic energy, is neg- 
ative under the barrier. In the classical manner this can 
be written in the form (compare with Eq. (jlip ) 



k r -\- k ? -\- 



1 



2y/r 2 + z 2 



1. 



(59) 



In our case the z axis is the tunneling direction and tun- 
neling probability is determined by Im k z . The smaller 



is this value the larger is the probability. Eq. ([59)) for- 
mally says that a reduction of Im k z is connected with 
increasing of Im k r . 

A localization near the z axis (the wave function expo- 
nentially decaying away from the axis) is equivalent to a 
large Im k r at that region and therefore leads to increase 
of tunneling rate along the z direction. A localization 
of the density around the z axis is possible but, as a 
"payment" for this, the wave function has a non-physical 
singularity on that line. Therefore one has to add the 
second branch to compensate the singularity. But the 
second branch exponentially increases away from the z 
axis which violates the construction. 

The above statement holds when some particular value 
of k r is chosen. If we consider a continuous set of waves 
with various k r , properties of the resulting wave function 
can be completely different. It is exactly our case. Due to 
an interference in the set of second branches, they com- 
pensate each other at finite z excepting a narrow (non- 
semiclassical) vicinity of the z axis. This occurs because 
of lack of a saddle point in the /3-integration. As a result, 
the state localized near the z axis (channel) is formed. 

The channel formation is strongly connected with cusp 
caustics. Without the Coulomb center the flux is shown 
in Fig. 03a). The Coulomb forces separate the flux by 
one, inside the caustic A%, and another, outside the caus- 
tic B\, as illustrated in Fig. Etb) . The former reaches the 
center but the latter passes far from it. Curves of con- 
stant \tp\ are almost conformal with the caustic A\ when 
they are close to it as in the inset in Fig. |4j The flux gets 
away from the cusp of the caustic A\ and proceeds to the 
center as shown in the inset. In this region is larger on 
the curves which are closer to the z axis. Therefore the 
value of on the caustic A\ is carried to the Coulomb 
center along the narrow channel containing the z axis. 
Not too close to the center the radius of the channel is 
proportional to l/B which is the de Broglie wave length 
in physical units. 

This is cusp driven tunneling when the cusp of the 
caustic Ai forms the non-small flux to the center. In the 
vicinity of the Coulomb center the flux passes the bottle 
neck of the size of l/B 2 which is Bohr radius in physical 
units. Then at negative z the flux proceeds to the left. 

We see that the cusp shape of caustics plays a crucial 
role. To get a cusp caustic one should have the certain 
form of a wave function at large distances. In our case the 
cusp position on the z axis, e 2 /(Esm 2 9o), is determined 
by the characteristic angle 9q of the incident conical flux. 
When the incident flux is just a plane wave, as in the 
conventional scattering problem, there is no cusp caustic. 
In this case only the usual caustic exist as in Fig. [T] 

To summarize the discussion of the phenomenon, one 
can propose its general qualitative explanation based on 
global properties of the flux. The Coulomb field separates 
the incident conical flux, with the characteristic angle 9q, 
according to the geometrical rule. Namely, all rays with 
larger angles are turned upward and ones with smaller 
angles are turned downward. This is equivalent to a sep- 
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aration of the two streams by caustics. In continuation 
toward the Coulomb center, the internal caustic should 
shrink to a point on the z axis before the center. Oth- 
erwise the flux classically reaches the Coulomb center. 
The caustic surface has to be of a cusp shape near this 
point to allow the flux to proceed from the caustic. An 
obtuse shape would result in a non-physical flux to this 
point on the surface which is locally perpendicular to the 
z axis. So the conclusion about the cusp formation fol- 
lows from general arguments without calculations. The 
essential condition is the conical incident flux but not a 
plane wave. 

The fusion probability is proportional to a neutron 
yield. The neutron yield of the process depends on the 
weight function cj(/3) (fl7|) and, therefore, on details of 
the incident flux. With the choice (1201 the yield is not 
exponentially small as exp(— 6174.9) (see Sec. [TTJ) but it 
can be "normally" small due to a preexponential factor 
in the tunneling probability. One can increase this factor 
by variation of a pre- factor in the weight function ([2D")) . 

We consider in the text the axially symmetric wave 
function of two deuterons related to the magnetic quan- 
tum number m — 0. A finite m reduces the effect due to 
formation of a centrifugal repulsive barrier. 

In this paper we briefly mention experimental schems 
for formation of the proper flux of particles resulting in 
the cusp phenomenon. 

One of experimental ways to produce the conical flux 
of deuterons ([56]) , shown in Fig. Ufa), is to confine them 
in a long tube with rigid walls, for example, in a nan- 
otube. Another way is to push deuterons (atoms) to pass 
through a diffraction grid of a conical shape. Since the 
de Broglie wave length is of the order of 1 A, one can 
use a natural crystal lattice. A setup with slits also can 
be used. This is a situation of quantum lens. We will 
discuss the details elsewhere. 

Above two bare deuterons are considered. From the 



practical standpoint it can be more convenient to deal 
with a substance (heavy water, for example) consisted 
of molecules with deuterium. Inside a single molecule of 
heavy water the deuterons are in the well with vibrational 
energy levels. An external laser radiation can influence 
a quantum state of deuterons in the well. One can put 
a question whether the radiation is able to create some- 
thing like cusp state in the well and which pulse shape 
should be used for this. 

Anyway a formation of a required particle flux in ex- 
periments needs a detailed study which is outside this 
paper. 



XI. CONCLUSIONS 

Fusion of two bare (supposed to be isolated from ev- 
erything) deuterons of room temperature energy can be 
possible. The penetration across the Coulomb barrier, 
which is the main obstacle for low-energy fusion, strongly 
depends on a form of the incident flux on the center at 
large distances from it. In contrast to the usual scatter- 
ing, the incident wave is not a single plane wave but a 
superposition of plane waves of the same energy and var- 
ious directions. As a result of interference of them, the 
wave function close to the Coulomb center is determined 
by the cusp caustic which is probed by de Broglie waves. 
The de Broglie waves get away from the caustic cusp and 
proceed to the Coulomb center providing the not small 
probability of fusion (cusp driven tunneling). 
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